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Conductance fluctuations in metallic nanogaps made by electromigration 
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We report on low temperature conductance measurements of gold nanogaps fabricated by controlled electro- 
migration. Fluctuations of the conductance due to quantum interferences and depending both on bias voltage 
and magnetic field are observed. By analyzing the voltage and magnetoconductance correlation functions we 
determine the type of electron trajectories generating the observed quantum interferences and the effective 
characteristic time of phase coherence in our device. 



I. INTRODUCTION 

Metallic nanogaps made by electromigration are key 
components of three terminal single molecule electronic 
devices. Besides their potential use as nanometer-spaced 
contact electrodes in future molecular devices, they are 
also ideal systems to investigate quantum transport prop- 
erties at the nanoscale. As in other mesoscopic systems, 
quantum interferences are thus expected to have a strong 
influence on electronic transport properties of nanogaps 
at low temperature. The typical signatures of these phe- 
nomena in diffusive wire are weak localization and uni- 
versal conductance fluctuations^— . Voltage dependent 
conductance fluctuations have been reported so far in 
related systems such as wires interrupted by a tunnel 
junction^ — or an atomic contacts^—, however they have 
not been investigated in metallic nanogaps. This is an im- 
portant issue for single molecule based transistors since 
the measurement and the operation of this kind of molec- 
ular devices can be hindered by quantum interference ef- 
fects taking place in the wires, which produces variations 
of the current unrelated to the specific signatures of con- 
duction processes through the single molecule. 

In this article we report conductance fluctuations in 
gold electromigrated nanogaps as a function of voltage 
and magnetic field. By analyzing the correlation func- 
tions of the fluctuations we determine the physical mech- 
anism responsible of quantum interferences and we ex- 
tract the effective characteristic time of phase coherence 
in our device. 



II. EXPERIMENT 

Samples are fabricated at room temperature in high 
vacuum (10~ 6 mbar) by controlled electromigration of a 
gold nanowire 500 nm long, 100 nm wide and 20 nm 
thick. The electromigration procedure has been pre- 
viously describe d 11 ' 12 , however the final steps leading 
to the nanogap formation have been changed following 
O'Neill et &\M. We initially thin the nanowire until the 
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FIG. 1. (Color online) Scanning electron microscope images 
of sample 1 (a) and 2 (b) after electromigration and thermal 
cycling. 



conductance reaches typical values of 3-4 Go, Go = 2e 2 /h 
being the quantum of conductance for degenerate spins. 
Then we take advantage of the high diffusion constant 
of gold atoms at room temperature^ by letting the 
nanowire spontaneously evolve while monitoring its con- 
ductance value. Atoms rearrangement is accompanied by 
conductance steps, indicating the evolution of the con- 
tact through a small number of conduction channels. In 
several hours a typical conductance value of 0.1 Go is ob- 
tained corresponding to bare nanogaps in the tunneling 
regime. In these conditions we cool down the sample to 
our dilution refrigerator base temperature (60 mK). The 
cooling down procedure typically decreases the conduc- 
tance by one order of magnitude to 0.01-0.05 Go- 

We present results on two samples obtained with the 
same procedure but with different final conductance 
mean values, 0.0086 Go for sample 1 and 0.050 Go for 
sample 2. Scanning electron microscope (SEM) images 
of the two samples are shown on Fig. 1 (a) and (b) 
respectively. Starting from the 500 nm long nanowire, 
electromigration induces an asymmetric breaking. We 
note that for sample 1 the resulting electrodes lengths 
are ~100 nm and ^400 nm (Fig. 1 (a)), and for sample 
2 we measure similarly ^200 nm and ~300 nm (Fig. 1 
(b)). Such asymmetry in electrodes length plays a funda- 
mental role in electron coherence as we will show in this 
article. Near-identical results for the conductance fluc- 
tuations as a function of magnetic field and bias voltage 
were also found on a third broken nanowire fabricated 
fr on the same substrate, showing the same assymetry and 
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having a conductance 0.0005 Go- Note that SEM images 
are taken after warming up the samples, as a consequence 
the nanogaps sizes visible on Fig. 1 are enlarged due to 
thermal cycling. 

Conductance measurements are performed using a 
lock-in amplifier technique with an AC voltage V m — 
30/j.V at a frequency of 77 Hz. Electrical leads are fil- 
tered by an RC distributed line thermally anchored to 
each stage of the dilution unit, connecting the dilution 
refrigerator mixing chamber to room temperature. In 
particular aim long RC distributed filter is located at 
the base temperature stage where the sample is mounted 
in a microwave-tight copper box. The sample is posi- 
tioned at the center of a superconducting magnetic coil, 
which allows to apply a perpendicular magnetic field. 



III. RESULTS 

Transport measurements on the examined samples 
show no signatures of Coulomb blockade or Kondo 
physics due to residual gold clusters ensuring that 
nanogaps are in a well established tunneling regime. 
When measuring the samples the voltage drop is localized 
at the nanogap, since its conductance is much smaller 
than the conductance of the leads (^75 Go). This allows 
to vary the energy interval over which electronic interfer- 
ences occur resulting in voltage dependent conductance 
fluctuations. 

Figs. 2 (a) and (b) show conductance measurements 
for sample 1 and 2 as a function of the bias voltage for 
different values of the magnetic field, from to 1.3 T. 
A non linear background due to the voltage dependence 
of tunnel transport has been subtracted. Conductance 
fluctuations are clearly evident and reproducible. Note 
the presence of a well defined dip at V = for all val- 
ues of the magnetic field. This zero-bias anomaly (ZBA), 
independent of the sample, can be explained by dynami- 
cal Coulomb blockade due to the dissipative electromag- 
netic environment surrounding the nanogap. In the in- 
set of Fig. 2 (b) we analyse quantitatively this anomaly 
for the G(V) data normalized to the tunnel conductance 
Gt at B = 1.1 T. The dissipative environment of the 
nanogap is modeled as a parallel capacitance and a sc- 
ries resistance^. For an uncontrolled environment, the 
scries resistance is expec ted to be close to the vacuum 
impedance Zo = \J Ho/eq at microwave frequencies. The 
fit, shown as solid line, is then performed using two free 
parameters, the junction capacitance G and the elec- 
tronic temperature T e , and fixing the series resistance 
R = Z =377 fl. The resulting value C = 10 ± 6 fF 
can be related by geometrical arguments to the capaci- 
tance towards the gate electrode. The found electronic 
temperature T e = 500 ± 100 mK is consistent with the 
temperature corresponding to the AC excitation voltage 
~ 360 mK, such strong excitation was needed to get ac- 
curate measurements. The uncertainties are quite high 
due to the difficulty to infer Gt accurately because of the 



12 



11 



; 10 




(a) 




(b) 



1.00^ 
g- 0.9 




65 



60 i 



55 



50 



4-4 
V(mV) 



4 
V(mV) 



FIG. 2. (Color online) Conductance of sample 1 (a) and 2 (b) 
as a function of the bias voltage in the range -4 mV < V <6 
mV for different values of the magnetic field (from bottom to 
top B=0, 0.1, 0.3, 0.5, 0.7, 0.9, 1.1, 1.3 T). Curves are verti- 
cally shifted for clarity. Inset: Dynamical Coulomb blockade 
filis* of the zero bias anomaly of the G(V) curve at B= LIT. 
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FIG. 3. (Color online) Conductance of sample 1 (a) and 2 (b) 
as a function of the magnetic field in the range 0.5T < B <5.5 
T for different values of the bias voltage (from bottom to top, 
V = 1.1, 2.7, 4.1, 6.1, 7.1, 8.1, 10.1 mV for sample 1 and V = 
1.1, 3.1, 5.1, 7.1, 9.1, 11.1 for sample 2. Curves are vertically 
shifted for clarity. 



conductance fluctuations. 

When measuring the conductance as a function of the 
magnetic field at different bias voltages analogous fluctu- 
ations occur. Figs. 3 (a) and (b) show the G(B) curves 
for different values of the bias voltage, from 1 to 11 mV, 
for sample 1 and 2 respectively. A linear background, 
where present, has been subtracted. Note that we have 
verified that the G(B) curves are symmetric with respect 
to B = 0. 
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IV. DISCUSSION 

The observed fluctuations are signatures of interfer- 
ence of electronic trajectories in metallic wires, which 
highlights the wave nature of electrons^. We analyze our 
data by calculating the correlation function normalized to 
the variance defined as F(AB, AV) = {SG(V, B)8G{V + 
AV,B + AB))/(6G(V,B) 2 ). In Fig. 4 (a) and (b) 
F(AB = 0, AV) is plotted as function of AV" for the 
two samples (ZBA has not been taken into account in 
data computation). The experimental data (dots) are 
obtained from measurements reported in Fig. 2 (a) and 
(b), consisting of a set of 1000 data points. Following 
the ergodic hypothesis, changing the magnetic field cor- 
responds to changing the impurity configurations in the 
wire, so data are averaged on all the explored magnetic 
fields values. The inset of Fig. 4 (a) and (b) shows 
F(AB, AV = 0) as function of AS for the two sam- 
ples at V = 1.1 mV, obtained from a set of 2500 data 
points in the interval 0.5<B<3 T. 

These data are compared with two theoretical mod- 
els, previously used for describing fluctuations in metallic 
wires interrupted by a tunnel junction^ and in atomic-size 
point contacts 9 . On the one hand, our device may behave 
as an atomic contact with low transmission coefficient 
since the spontaneous breaking of the nanowire occurs 
via discrete conductance steps indicating a small num- 
ber of conduction channels. On the other hand, during 
the cool down procedure we measure a further conduc- 



tance decrease without characteristics steps. This reveals 
atomic rearrangement which may lead to a large number 
of low transmitted conduction channels, a configuration 
better described by the tunnel junction model. 

In these models the energy of the injected electrons 
can be tuned by the applied bias voltage. Impurities 
in the leads allow diffusion of electrons through multi- 
ple scattering processes, so that once an electron is in- 
jected in the lead, it can be backscattered. In the tun- 
nel junction model, conductance fluctuations are due to 
closed loop trajectories involving the two leads and in- 
cluding the junction. In the case of an atomic contact, 
conductance fluctuations are generated by interferences 
between electron waves directly transmitted by the con- 
tact and the waves backscattered by the diffusive leads 
and reflected at the contact towards the same lead. In 
both cases interference loops have in first approximation 
typical area equal to L x W, where W is the width of 
the electrodes and L is the characteristic length limiting 
electron phase coherence. A change AB in the magnetic 
field applied perpendicularly to the sample adds a phase 
shift Aip = 2ttABLW/(/)o to the electron wavefunction, 
where (f>o = h/e is the flux quantum. The typical field 
scale of fluctuations is given by 4>o/(LW). In the same 
way, a phase shift can be induced by an energy change 
eAV giving Aip — eAVL 2 / (hD), where D is the diffusion 
constant of electrons in gold. The corresponding voltage 
scale of fluctuations is hD/(eL 2 ). The full theoretical 
expression for the correlation function in the two cases^ 
is given respectively by 



(SG(V, B)6G(V + AV, B + AB))\ = J.^J / P 2 cos(^AW) e - 2t ^ e -*/ T - J 2 {\V m t) dt (1) 

(5G(V,B)5G(V + AV,B + AB)) 2 = (^) 2 ^±_Jl f° (-) 2 P d>2 (t) cos^AF^e-^e"'/™ J? ( \v m t) dt (2) 

h V£ J t h h 



where the index 1 and 2 refer to the tunnel junction and 
the atomic contact model. Here p is the electrode den- 
sity of states, T the atomic contact transmission coeffi- 
cient, V the DC bias voltage, tb the dephasing time due 

to the magnetic field, given by tb = 12— ^^r-, and r 
is the time limiting electron coherence. The first order 
Bessel function Ji(i V m t) takes into account the dephas- 
ing due to the amplitude of the applied AC voltage V m 
which has to be considered if eV m > Ujr. The function 
Pci,i(t) is the classical probability for an electron tun- 
neling through the nanogap to be backscattered towards 
the nanogap again for the two models (i = 1,2). In the 
geometry and energy range considered here, transport is 
quasi one-dimensional thus P c i,i(t) — J^ i Dt > where Ki is 
a factor including constant terms and geometrical details. 
The return probabilities appear with different power laws 



in Eq. (1) and (2). This is related to the particular 
electronic trajectories generating quantum interferences 
considered in the two models. In the tunnel junction 
model, the transport is supposed to be diffusive over the 
whole length of the wire except at the tunnel junction. 
Therefore, the diffusive trajectories contributing to con- 
ductance fluctuations comparable with the junction con- 
ductance, are closed loops crossing the tunnel barrier and 
involving both the left and right electrodes. As a conse- 
quence, assuming the diffusion constants on both sides 
are equal, Eq. (1) presents a square dependence on the 
classical return probability P c i,i(t). In the atomic con- 
tact model, the system is described as a central ballistic 
region containing the atomic contact surrounded by dif- 
fusive regions on either side. The linear dependence on 
the return probability in Eq. (2) can be obtained by 
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FIG. 4. (Color online) Normalized correlation function (blue 
dots) for sample 1 (a) and 2 (b) as a function of AV. Insets: 
normalized correlation functions as a function of AB. Data 
are compared with the tunnel junction model (black curve) 
and the atomic contact model (red dashed curve). 



the length over wh ich e lectron coherence is preserved by 
the formula L = yDr, giving L ~ 110 nm. This result 
is consistent with the characteristic length extracted by 
estimating AV from data in Fig. 2 (a) and (b) using 
the formula = hD/{eL 2 ). We obtain L « 100-140 
nm with AV » 1-2 mV. These values are lower than the 
expected phase coherence time and length in gold, 
of the order of 200-300 ps and 20-30 /im respectively at 
the experimental temperature. This is not surprising by 
considering the typical dimensions of our devices. Elec- 
tron coherence in fact is limited by the length of the 
nanogap electrodes and the large reservoirs. As shown in 
Fig 1, initially the nanowire is 500 nm long, after electro- 
migration the nanogap is formed not in the center of the 
wire, the shorter electrode has length as small as 100 nm. 
Note also that L remains smaller than the thermal length 
L T = y/hD/k B T e = 280 nm for T e = 500 mK, indicating 
that electron coherence is not temperature limited. 



V. CONCLUSION 



considering electron waves being transmitted by the con- 
tact and injected from the central ballistic part in one 
of the diffusive regions. They can be backscattered with 
probability P c i,2(t), reflected by the contact, and inter- 
fere with the directly transmitted waves, thus generating 
conductance fluctuations. 

Full and dashed lines in Fig. 4 represent the fit to the 
experimental data on the basis of the two models. When 
fitting the voltage correlation function we set W = 10~ 7 
m and S — 2xl0~ 15 m 2 according to the geometry of 
the samples. Moreover, since, at AB = 0, the calculated 
voltage correlation functions weakly depend on the dif- 
fusion constant, we take D = 0.005 m 2 /s estimated from 
the Einstein relation^ by measuring the nanowire resistiv- 
ity before electromigration at room temperature^. The 
only free parameter is then the time r. Clearly the tunnel 
junction model seems more appropriate to describe our 
experimental data. This result indicates that even if dur- 
ing the fabrication process, the nanogap evolves through 
a small number of conduction channels, final atoms re- 
arrangement at low temperature corresponds better to 
the tunnel junction limit of a large number of conduc- 
tion channels with low transmission. Moreover the tun- 
nel junction model suggests that interferences in our de- 
vices are due to trajectories involving the two electrodes 
and the nanogap. Using Eq. (1) we extract r obtain- 
ing similar results for the two samples, r ~ 2.3±0.3 ps. 
Unlike the voltage correlation functions, the calculation 
of the magnetoconductance correlation functions at AV 
= strongly depends on D through the exponential de- 
pendence in tb- In order to fit the experimental data 
(inset of Fig. 4 (a) and (b)), we fix r at the value ex- 
tracted from the former fit and we leave D as free pa- 
rameter. We obtain D = 0.006±0.002 m 2 /s in agreement 
with our previous estimate. Knowing r we can calculate 



In conclusion, we have measured conductance fluctua- 
tions as a function of the bias voltage and the magnetic 
field in gold electromigrated nanogaps at low temper- 
ature. By analyzing the correlation functions on the 
basis of existing theoretical models, we conclude that 
the nanogap can not be described as a low transmit- 
ted atomic contact but rather as a tunnel element inter- 
rupting a diffusive wire. Fluctuations are due to quan- 
tum interferences related to closed electronic trajectories 
limited by sample geometry and including the nanogap. 
Since, in our device, electron coherence is limited by the 
length of the nanogap electrodes, increasing the length of 
the wire used to form the nanogap would increase the ef- 
fective characteristic length of phase coherence and even- 
tually would allow to probe the intrinsic electronic coher- 
ence length. This study shows that electron coherence in 
contact electrodes has to be taken into account in future 
single molecule electronic devices. Especially the effect 
of unwanted conductance fluctuations may be minimized 
by using short and wide source/drain electrodes in the 
design of three terminal single molecule devices. 
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